Physical realization of Photonic Klein Tunneling 
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General physical conditions for the occurrence of photonic Klein tunneling are studied, where (con- 
trolled) spontaneous emission from the devices considered plays a key role. The specific example 
of a simple dielectric slab bounded by two dielectric half spaces with arbitrary refractive indices is 
worked out quite in detail, the measured reflection and transmission probabilities being calculated 
analytically. It is found that, in given cases, the measured reflection probability may be arbitrarily 
large (for large incident wavelengths) irrespective of the fact that the transmission probability is 
exponentially suppressed or not. Other interesting features of photonic Klein tunneling driven by 
(controlled) spontaneous emission are as well envisaged for practical applications. 
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The experimental observation [l[ of the intriguing prop- 
erties of graphene Q, a genuine two-dimensional mate- 
rial composed of carbon atoms forming a honeycomb lat- 
tice, has recently drawn a very vivid interest by many 
physicists (see [1] and references therein). The main 
reason relies in the fact that the forbidden energy gap, 
through which the charge carriers tunnel, is suppressed 
to zero in graphene, differently from other materials, 
and, in addition, those peculiar charge carriers behave 
as Dirac fermions Q, Q. This means that electrons 
in such unique material behave effectively as massless 
Dirac fermions, thus allowing the possible experimental 
study of a number of key phenomena Indeed, the 
scattering of relativistic fermions, being described by the 
Dirac equation, is fundamentally different from that of 
non-relativistic ones. For example, relativistic electrons 
scattered off a potential step, with height higher than 
their energy, at normal incidence may exhibit a non-zero 
transmission probability, this being in sharp contrast to 
the intuitive result of non-relativistic quantum mechan- 
ics. Such aphenomenon, in general referred to as Klein 
tunneling |fj(, that is the below-barrier particle tunnel- 
ing without the exponential damping expected for non- 
relativistic particles, has never been observed experimen- 
tally, since too much strong electric fields are required to 
observe it in elementary particles. Nevertheless, Klein 
tunneling may be realized on graphene sheets Q: since 
a sufficiently strong potential, though repulsive for elec- 
trons, is attractive for holes, it gives rise to hole states 
in the barrier to form channels through which electrons 
can penetrate the barrier Q. Several experiments have 
been made to observe such phenomenon, and we refer the 
interested reader to the existing literature [7] for further 
discussion. 

Very recently, even in view of possible intriguing appli- 
cations, optical analogues of Klein tunneling have been 
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as well proposed [8j-[13j, thus capturing the interest of 
many. This is the case, for example, of light propagation 
in deformed honeycomb photonic lattices [H [l(|, whose 
band structure is similar to that of graphene, or light 
refraction at the interfaces between positive and nega- 
tive index materials fill ]. Other proposals regard spatial 
light propagation in binary waveguide arrays [lOj, one- 
dimensional stationary light pulses in an atomic ensemble 
with electromagnetic induced transparency (in the limit 
of tight spatial confinement) 1121 . or Klein tunneling of 
light in fiber Bragg gratings [13| . 

In these papers, the investigations on photonic ana- 
logues of relativistic tunneling phenomena are carried out 
exactly as made for an alog ues of non-relativistic tunnel- 
ing (see, for example, [14| and references therein), that 
is by writing down the dynamical equations for the light 
propagating in the device considered and then exploiting 
the formal mathematical analogy between such equations 
and the Dirac equation for electrons. For example, in 
]13j it is shown that the counter-propagating waves in a 
fiber Bragg grating satisfy coupled-mode equations that 
can be exactly cast in the form of a Dirac equation in 
presence of an electrostatic field, provided a suitable def- 
inition for a "spinor" wavefunction. The optical analogue 
of the forbidden energy region is given by the photonic 
stop band of the periodic grating. Pulse propagation in a 
fiber Bragg grating with a suitably designed chirp profile 
can then be used to mimic the relativistic tunneling of a 
wavepacket in a potential step. Other, different realiza- 
tions, as mentioned above, are treated analogously. 

Such a way of reasoning is certainly legitimate, and 
leads to correct results. Nevertheless, by focusing mainly 
on the formal analogy among mathematical equations 
does not provide a full physical insight in the true phe- 
nomenon considered, and thus only particular cases are 
usually studied. In this note, instead, we will discuss 
some general conditions for the realization of photonic 
Klein tunneling, mainly focusing on the physical inter- 
pretation of it. 

Let us start with the typical example of an electron 



FIG. 1: A diagrammatic representation of Klein tunneling of 
electrons by an electrostatic potential source. 



source 
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(with relativistic energy) impinging on a potential bar- 
rier of height Vq created, for example, by a sufficiently 
strong electric field. As well known, Klein tunneling takes 
place only when Vq is larger than twice the particle rest 
energy 2m e c 2 , so that pair production can occur, accord- 
ing to the diagrammatic representation [26] in Fig. [I] 
For photonic tunneling, such a limitation is not present, 
since photons have no rest energy, so that, in principle, 
photonic Klein tunneling should occur more easily. This 
simple expectation, however, confronts with the effective 
realization of the potential barrier for the photons or, in 
other words, with the coupling of the incoming photon 
to the barrier: for electron tunneling, indeed, electrons 
couple directly to the electric field originating the poten- 
tial barrier. A naive transposition of electron tunneling 
to photonic tunneling would proceed according to the 
diagrammatic representation in Fig. [2^, where photon- 
photon coupling even in vacuum is indeed possible in a 
full quantum picture (at least at one loop level) through 
the creation of virtual electron-positron pairs, when non- 
linear effects come out. In such a case (photons impinging 
on an "electromagnetic barrier"), Klein tunneling effec- 
tively reduces to the phenomenon of photon splitting, ex- 
tensively considered in the literature 
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(and similarly for the B field), where the "background'' 
(source) quantities are given by 

P = 2C [2 (E 2 - c 2 B 2 s ) E ; 



M = -2c 2 C [2 (E 2 - c 2 B 2 



f 7c 2 (E s -B s )B s ] (2) 
B S + 7(E S -B S )E S ] (3) 



with C = 2a £qTi /4:5m^c 5 (a = e /2ehc is the fine struc- 
ture constant), under given conditions allow the photon 
propagation in a way similar to that in an undersized 
waveguide, thus mimicking a tunneling phenomenon. 
However, as largely debated in the literature 
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and as evident from the equations above, such effect is 
almost unobservable in laboratory experiments and cer- 
tainly not of practical use in appliction, due mainly to 
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FIG. 2: A diagrammatic representation of photonic Klein tun- 
neling: a) by an electromagnetic potential source; b) by an 
electron source (material medium). 



the extremely small value of the effective coupling con- 
stant C, while it may play a non-negligible role in several 
astrophysical environments. For this reason, we do not 
consider here such a possibility. 

An alternative comes from the coupling of the tun- 
neling photon to an "electron" (rather than electromag- 
netic) source, as diagrammatically depicted in Fig. [5b; 
that is, the potential barrier is given by electrons in a 
suitable material, such as a refractive medium. This is 
not an example strictly equivalent to electron tunneling, 
but is a common and easy way to achieve (even standard, 
not Klein) photonic tunneling. In such a case, referring 
to Fig. [2}d, the additional photon |27j is produced by the 
electrons in a given medium by means of spontaneous 
emission, which is, again, a pure quantum effect. The 
physical realization of this version of photonic Klein tun- 
neling may proceed through different mechanisms, but 
its practical achievement is related to the possibility of 
effective (and controlled) spontaneous emission. 

As an example, we may consider the spontaneous emis- 
sion of photons by impurity atoms or molecules in ho- 
mogeneous dielectrics. The vacuum emission rate of a 
randomly oriented dipole with moment d is given by fri\ 



d 2 ul 



(4) 



for a transition of frequency ljq, but, as it is now well 
established, the rate at which an excited species under- 
goes spontaneous emission depends on its environment 
[18| [19j . The spontaneous emission rate in a bulk di- 
electric, for instance, is scaled by the real part of the re- 
fractive index n = yfe of the medium at the frequency of 
the transition In general, the spontaneous emission 
field can be strongly modified if the emitting atom is sur- 
rounded by materials of different composition and shape, 
since the existence of material boundaries (of any kind) 
in the vicinity of the radiating species simply changes the 
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FIG. 3: Photonic Klein tunneling through a dielectric slab. 



strength and distribution of the electromagnetic modes 
with which the emitter interacts, then resulting in an al- 
tered spontaneous emission rate. In particular, when an 
increase in the density of states occurs, the rate of spon- 
taneous emission can be enhanced over the free space 
value (@|. Similar results are obtained when the atoms 
are placed between mirrors or in resonant cavities [l8j], 
in which case larger modifications are possible. 

Spontaneous emission rates are also altered in photonic 
crystals [2(| , where the dielectric constant is periodically 
modulated on length scales of the wavelength of light. 
Due to multiple Bragg reflections, light in a certain fre- 
quency range cannot propagate in the crystal, so that 
the radiative density of states in this photonic band gap 
is zero and, therefore, spontaneous emission of excited 
atoms embedded in these materials is expected to be in- 
hibited. However, at other frequencies it has been shown 
that spontaneous emission can be either enhanced or re- 
duced, so that in photonic band gap materials a nearly 
complete control over spontaneous emission rates might 
be achieved. This opens the intriguing possibility to have 
a single material where standard and Klein tunneling 
may be observed for incident photons with different fre- 
quencies. 

A number of other devices (left-handed materials [2l[ , 
sharp metallic tips |22j |. waveguides (23[, quantum dots 
24], etc.) may as well be used to modify the emission 
features, all of them driven by the zero-point fluctuations 
of the electromagnetic field. It is at last clear that the 
key ingredient relevant for the practical realization of the 
photonic Klein tunneling we are dealing with is just the 
general result that spontaneous emission of an emitter 
can be controlled to some extent, and then engineered by 
tailoring the surrounding structure on a transition wave- 
length scale. 

Let us now address the question of finding general con- 
ditions for photonic Klein tunneling to occur. In order 
to get definite results, we shall consider a specific ex- 
ample for illustrative purposes, that is a dielectric slab 
bounded by two dielectric half spaces with arbitrary re- 
fractive indices, as shown in Fig. [3J with the incident 
wave coming from the negative x region. The photonic 
propagation may be described simply by a scalar field, 



which we denote with ip in analogy to the Dirac case: 
it corresponds to some scalar component of the electric 
or magnetic field satisfying the appropriate electromag- 
netic Hclmholtz equation. In the most general case, by 
allowing non-vanishing spontaneous emission in the two 
barrier-free regions (in both negative and positive x di- 
rections), the solution of such equation can be written 
as 

!V>i(x), for x < 0, 
i/)o(x), for0<x<a, (5) 
ip 2 (x), for x > a, 

with 

^i(ar) = [l + 5i(a;)] e^ x 
+ {R[1 + St(x)} + + TS 2 (x)} e-^ x , 

Mx)=Ac- x " x + Be X0X , (6) 

M*) ^S 2 (x)e- tk ^ x -^ 
+ {T [I + Sx{x)] + S 2 (x) + RS 2 {x)} e lk ^ x - a \ 

where k\ , k 2 are the real wavevectors in the barrier- free 
regions, i\o is the imaginary wavevector in the < x < a 
region, and i?, T are the reflection and transmission co- 
efficient of the barrier, respectively. The real source 
functions Si(x), S 2 {x) take into account the spontaneous 
emission rates (and, if any, possible absorptive effects): 
for a dielectric slab as that considered here, these have 
been calculated in [25j as functions of the position, so 
that we refer the reader to this paper for further dis- 
cussion. However, since we are interested in general 
conditions, the explicit (though numerical) expressions 
for these are not really necessary here, as we shall see 
below. Note that, in tpi(x), the 5i(x)-terms accounts 
for the (isotropic) spontaneous emission in region 1, the 
R Si (x)-term accounts for the reflection of such signal 
by the barrier, while the T S 2 (x)-term accounts for the 
spontaneous emission signal in region 2 transmitted by 
the barrier in region 1. Analogously, in ^ 2 {x), the S 2 (x)- 
terms accounts for the (isotropic) spontaneous emission 
in region 2, the R S 2 (x)-teim accounts for the reflection 
of such by the barrier, while the TSi(a;)-term accounts 
for the spontaneous emission signal in region 1 transmit- 
ted by the barrier in region 2. 

The matching conditions at interfaces are given by 

M0) = M0), ^(o) = Vo(o), (7) 

00 (a) = "02(a), ip' (a) = i>' 2 (a) , 

where the prime denotes differentiation with respect to x. 
By substitution into Eq. (J6j> we can find the explicit ex- 
pressions for the function and, then, for the reflection 
and transmission probabilities. However, the resulting 
expression for such general case are excessively lengthy 
and no useful physical insight can be gained easily. We 
instead consider few interesting special cases, relevant 
for practical applications, the physical meaning of the 
corresponding results being particularly transparent. In 
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addition, we consider a symmetric barrier, with k 2 = ki 
(for k 2 k\ the calculations are a bit more involved, 
while no further physics comes out), in the fully opaque 
limit xo a ~~ ^ 00 ■ The last assumption is justified when 
focusing on the dominant terms of the Klein tunneling, 
without the characteristic exponentially decreasing fac- 
tors e -2x ° a . Of course, corrections including such factors 
are always present, but here we neglect their effect by 
considering only very opaque barriers. 
First case. Let us assume that spontaneous emission 
takes place only in region 1, and only along the positive 
x direction, as for the incident signal. In such a case, 
ipi{x) and ip 2 (x) in Eqs. ([6]) reduce to: 



[1 + Si(x)] e lklX +R[1 + Si(x)] e~ lklX , 
T[l + Si(a-x)]e ikl{x - a \ 



(8) 



We have, then, just a rescaled tp function in the barrier- 
free regions with respect to the standard case (with no 
spontaneous emission). As expected, by imposing the 
matching conditions ([7]), in the approximations above we 
have the following results for the reflection and transmis- 
sion probabilities by the barrier: 
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ITT 



0. 



(9) 



However, these probabilities are not directly obtained in 
experiments, since the measured probabilities are nor- 
malized with respect to only the signal emitted by the 
source employed and not to the total signal impinging on 
the barrier (including the spontaneous emission source). 
In the present case, then, the measured reflection and 
transmission probabilities are evaluated from the ratios 
| J R| 2 [1 + S'i(a;)] 2 and \T\ 2 [1 + Si (a - x)} 2 , respectively. 
Given the effective x dependence, these ratios should be 
integrated over the effective measurement regions but, 
taking into account the fact that the spontaneous emis- 
sion rate is dominant near the edges of the barrier [25[ , as 
a good approximation we can evaluate the source terms 
above just at the corresponding edges (28j, thus obtaining 
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IT, 



measured 



0. (10) 



From these results we get the quite obvious conclusion 
that, in the present case and in the approximations 
above, the reflection probability is greater than one due to 
the spontaneous emission, while the transmission proba- 
bility is exponentially suppressed. 

Second case. Let us now assume that spontaneous emis- 
sion takes place again only in region 1, but only along the 
negative x direction, as for the incident signal reflected 
by the barrier. In such a case, fa[x) and ip 2 (x) in Eqs. 
([6]) reduce to: 



fa{x) 



_ Akxx 



(/> 2 (i) = re*( I -°). 



+ [R + S 1 {x)]e~^ x , 



(11) 



In the same approximations as above, the measured re- 
flection and transmission probabilities are now given by 



|-R + Si(0)| 2 and \T\ 2 respectively, so that, by imposing 
the matching conditions ([7]), we get: 
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|2 k\ + [ X0 + suo)r 



measured 



measured 
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Again, then, signal transmission is exponentially sup- 
pressed, but now the reflection probability exhibits an 
interesting property: it is greater than one if S[ (0) is pos- 
itive, while it keeps lower than one in the opposite case. 
In other words, the super- unitary reflection is determined 
by the sign of the slope of the spontaneous emission rate 
which, as deduced from the results reported in Ref. [25j . 
depends on the characteristics of the apparatus employed 
(refractive indices of the slab, emission wavelength and 
thickness of the slab). 

Third case. The two particular cases just considered are 
generalized by assuming that the spontaneous emission 
(taking place only in region 1) is isotropic, that is it takes 
place along both the positive and negative x directions. 
Now, then, fa{x) and ip 2 {x) in Eqs. ([6]) reduce to: 

i>!{x) = [1 + Si(a;)]e lfcia: 

+ {R[l + S 1 {x)]+S 1 (x)}e- tk ^, (13) 
ifa(x) = T[l + Si (a - x)} e lkl(x - a l 

The measured reflection and transmission probabilities 
are given by \R [1 + Si(0)] + Si(0)| 2 and |T [1 + Si(0)]| 2 
respectively, and, by imposing the matching conditions 
(UJ), we have: 



\R 



k 1+ xo + SM 



measured 



—2 



*i +Xo 



[1 + Si(0)], (14) 



IT, 



measured 



0. 



where fei = fc x [l + Si(0)], x = Xo[l + Si(Q)] + S[(0), 
S[(0) = S[(0)/[l + Si(0)]. Even in the present case, 
then, the super-unitary reflection depends on the sign 
of the slope of the source term; in particular, again the 
reflection probability is certainly greater than one for 
S[(0)>0. 

Fourth case. As a final example, let us consider the case 
with the spontaneous emission taking place only in re- 
gion 2, along the negative x direction (of course, there is 
no interesting result for the reflection probability when 
light is spontaneously emitted only along the positive x 
direction, while changing the transmission properties in 
a fashion similar to that of the second case above). Now, 
tpi(x) and ip2(x) in Eqs. ([S]) reduce to: 



fa ( x ) = e lklX + [R + TS 2 (a- x)] e 
fa(x) = [T + RS 2 (x)]c tkl( - x - a '> + S 2 {x)e lk ^ x - a \ 



(15) 



The measured reflection and transmission probabilities 
are given by \R + TS 2 {a)\\ 2 and \T + RS 2 (a)\ 2 respec- 
tively. Differently from the previous cases, the trans- 
mission probability is now not exponentially suppressed; 
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indeed, by imposing the matching conditions ([7]), we get: 



I -^measured | 



IT, 



measured | 



[(fc?+^)-^ 2 ] 2 + 45 2 2 ( fl )^ 2 (a)fc? 

[{kl+xl)~^] 2 +AS'i{a)kl : 
Sj{a) m+xl)-K 2 ] 2 + 4S 2 2 («)fc 2 
[{kl+xl)-K 2 ] 2 +±S'i{a)kl 



(16) 

where k 2 = Sf(a)fc? + [S 2 (a)xo ~ S' 2 {a)] 2 . Note that 
both the reflection and the transmission probabilities are 
greater than one (as expected) for any value of S 2 (a) 
and value and sign of S' 2 (a). The total probability is, of 
course, again greater than one and, differently from the 
second and third case above, does not depend explicitly 
on /ci,xo> the following relation holding true: 

(-^measured) ~t~ |^mcasurcd| — 1 + 5*2 {a) • (17) 

A more general case, which is a better representation of 
physical reality, would include spontaneous emission from 
both barrier-free regions but, in such a case, mathemat- 
ical results become less transparent. Nevertheless, even 
in this general case the reflection probability is greater 
than one, while tunneling is not exponentially damped. 

From the results obtained above, it emerges quite clear 
that photonic Klein tunneling is not so much due to the 
mathematical structure of the electrodynamic equations, 
similar to the Dirac equation for electrons, but rather 
the physical content behind them. For electrons, this 
translates particle- antiparticle pair creation from the po- 
tential barrier. For photons, particle production from 
an electromagnetic potential barrier is as well a possi- 
ble phenomenon due to effective non-linearities in the 
Maxwell equations arising from quantum loop corrections 
(the phenomenon reduces, in a sense, to the well-known 
photon splitting); however, the extremely smallness of 



such corrections renders the effect of no practical use in 
applications. Instead, general Klein tunneling of pho- 
tons may be easily realized when they impinge on an 
"electron" (rather than electromagnetic) barrier, that is 
a refractive medium, and spontaneous photon emission 
takes place. This effect is of somewhat practical interest 
since spontaneous emission, contrary to naive expecta- 
tions, can be controlled in such media (for example, in 
photonic band gap materials) to some extent, as recalled 
above. Interestingly enough, as shown in the specific ex- 
ample of a dielectric slab considered here, photonic Klein 
tunneling may be realized even independently of the fea- 
tures of the spontaneous emission sources (fourth case). 
With the exception of the first case considered above, 
where spontaneous emission takes place only in the in- 
cident signal region along the direction of the incident 
signal, the measured reflection probability depends on ki 
and XO) and may be made arbitrarily large for large inci- 
dent wavelengths (and small imaginary wavevector of the 
barrier medium). Equally interesting, likely even in ap- 
plications, is the prediction that the measured reflection 
probability may be lower than that in the standard sce- 
nario (no Klein tunneling) if spontaneous emission takes 
place in the incident signal region along the direction op- 
posite to that of the incident signal (second and third 
cases), for given features of the source (negative slope of 
the spontaneous emission rate). 



The results obtained here for a simple dielectric slab 
may be generalized to different devices, so that a variety 
of possible implementations of photonic Klein tunneling 
may be realized in several applications, provided that 
the corresponding spontaneous emission properties can 
be easily managed. Such a challenge then claims for more 
theoretical and experimental studies, which will certainly 
lead to novel applicative solutions in the very near future. 
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